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ABSTRACT 



We present the results of a reverberation mapping (RM) campaign on the 
black hole (BH) associated with the active galactic nucleus (AGN) in SDSS 
J114008.71+030711.4 (hereafter GH08). This object is selected f rom a sample 



of 19 candidate intermediate mass BHs (M B h < 10 6 M Q ) found by lGreene fc Ho 



f l2004h in the Sloan Digital Sky Survey (SDSS). We used the Hobby-Eberly Tele- 
scope to obtain 30 spectra over a period of 178 days in an attempt to resolve 
the reverberation time lag (r) between the continuum source and the broad line 
region (BLR) in order to determine the radius of the BLR (Rblr) in GH08. We 
measure r to be 2 days with an upper limit of 6 days. We estimate the AGN 
luminosity at 5100 A to be AL 510 o ~ 1 1 x 10 43 erg s -1 after deconvolution from 
the host galaxy. The most well calibrated Rblr - L relation predicts a time lag 
which is 4 times larger than what we measure. Using the measured H/3 full- 
width-at-half-maximum of 703 ± 110 kms -1 and an upper limit for Rblr = 6 
light days, we find Mbh ^$ 5.8 x 10 5 M as an upper limit to the BH virial mass in 
GH08, which implies super— Eddington accretion. Based on our measured Mbh 
we propose that GH08 may be another candidate to add to the very short list of 
AGNs with M BH < 10 6 M determined using RM. 
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1. Introduction 



The determination of fundamental physical parameters is paramount for a complete 
understanding of the mechanism powering Active Galactic Nuclei (AGNs) and their role in 
galaxy evolution. The evidence for coevolution of the host galaxy and the central massive 
b lack hole is found observa tionally by the well known M— cr* relation (IFerrarese fc Merritt 



2000; Tremaine et al. 



2002) b ut it's applicability for low mass sources is still under debate 



(see for e.g., iKormendy et al.ll201ll ; iMathur et al.ll201ll ). AGN come in a variety of flavors 
and are found to have a range in bolometric luminosity from as low as 10 42 erg s -1 to beyond 
10 48 erg s _1 , and accordingly, have a large range in their central black hole mass (Mbh) from 
around 10 5 M up to several 10 9 M . For about 20 years the method which has yielded 
the most robust determinations of Mbh in AGN, over a range in luminosity and Mbh, is 
reverberation mapping (RM). The RM technique allows determination of the time-lag (r) 
between a change in continuum emission from an accretion disk and the delayed response due 
to light travel time of line emission from high velocity clouds in a distant broad-line region 
(BLR). Combining measurements of the velocity range {AV) from an emission line width 
originating in the BLR and the distance of the BLR from the BH (Rblr = cr, where c is the 
speed of light) gives the virial mass within the BLR. Since the BH dominates the mass inside 
the BLR, M B h = /AX^RblrG" 1 , where / is a scaling factor that depends on the geometry, 
kin ematics and orientati on of the BLR and G is the gravitational constant (for a full review 



sec 



Peterson et al.l 12004 ). Assuming an isotropic BLR circular velocity field and that the 



full w idth at half maxim um (FWHM) reflects two times the typical BLR velocity yields / = 
0.75 (IKaspi et al.l 120001 ). The most appropriate value of / to use depends on the method 
used to find AV and reported values vary from . 5 up to 5.5 depending on the FWHM of 



the line and the FWHM/<7zj ne ratio. ICollin et al.l ( 120061 ) define 2 populations based on the 



FWHM/<7/j ne ratio being below and above 2.35 (the intrinsic ratio when fitting lines with 
a Gaussian) and estimate / for each population to be 1.5 and 0.5 respectively. Therefore 
a reasonable / factor to use when measuring a line width measured from the FWHM of a 
Gaussian fit is 0.75. When Rblr is in light days and AV measured by the FWHM in 1000 
kms -1 , the BH mass in solar masses is: 



M BH = 1-464 x lO 5 i? BL RAK 2 M (1) 

When the RM determined Rblr is plotted against the AGN continuum luminosity a 
general trend is revealed where Rblr oc L a , the so called radius— luminosity (Rblr— L) power- 
law relation. An important caveat of the Rblr - L relation is the difficulty to accurately 
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deter mine the contributi on from the inner region of the host galaxy to the luminosity of the 
AGN. lBentz et al.l ( 120061 ) show, using high resolution Hubble Space Telescope (HST) images, 
that an accurate determination of the AGN luminosity with respect to the host galaxy 
luminosity can change the slope (a) of this relation significantly. A physical consequence of 
this relation is that for AGNs with very high luminosities, r can be very large, on the order 
of several to tens of years. Since the RM technique requires a high degree of accuracy in the 
calibration of spectra over an extended period of time for the most luminous AGNs, only a 
few of the highest luminosit y, and therefore the highest mass AGNs have been attempted 
with RM (IKaspi et al.l 120071 ). On the other hand, for AGNs with the lowest luminosity, r 
can be quite short, on the order of a few hours to a few days. Again, in terms of observations 
RM turns out to be difficult. If r is on the order of hours, then semi- weekly observations 
will not resolve the time lag. The largest group of AGNs with M B h determined using R M is 
the nearby Seyfert galaxies ( Wandel et al.l Il999t iBentz et al.l |2010| ; iDenney et al.l |2010| ) and 



several PG quasars kaspi et alil2000h who have r on the order of a few to several hundreds 
of days and 10 6 < M BH /M < 10 9 . AGNs with M BH < 10 6 M Q will have small r, on the 
order of hours to days, and slower bulk motion of clouds in the BLR as determined by the 
width of the broad lines. For example if r w 7 days then AV < 1000 kms -1 for Mbh < 10 6 
M Q . Only a very few AGN meeting these criteria have had Mbh robustly determined using 
R M The AGN in NGC 4 395 is one such object whose M BH ~ 3.6 x 10 5 M was determined 
by IPeterson et al ( 2005 ) using HST UV spectra to find a Civ time lag of about 1 hour. 
Desroches et al.l (120061 ) tenuously measure a ground based Balmer line time lag for NGC 
4395 that is consistent with the HST measurement. NGC 4395 is notorious in that it does 
not have a classic bulge and the refore does not follow the M— a* relation. NGC 4051 was 
found by Ishemmer et aD (12003^ to host a NLSl-like AGN with M BH « 5j| x 10 5 M Q 



measured using RM. IWandell (120041 ) subsequently show that NGC 4051 follows the M— cr* 
relation but h as an underluminou s bulge when compared to more standard type 1 AGN. 
More recently E)ennev et all (bold ) find M BH = (1.7 ± 0.5) x 10 6 M for NGC 4051 raising 
the initial estimate but still within the low mass regime. Low mass AGNs like NGC 4395 
and NGC 4051 are useful for constraining evolut ionary models since th ey may challenge in 
different ways the current coevolution paradigm (IKormendy et al.ll201ll ). 



The robust RM measurements from the moderate luminosity AGNs and the scaling 
relations derived from it, like the R B lr— L relation, can be directly applied to single spectra 
to determine M B h for large samples of AGNs. For this so called 'single epoch' method, 
M BB can be calculated using only 2 measurements, the AV of an emission line and the 
continuum luminosity. For lower redshift samples AV^ is usually found using the H/3 line and 
the luminosity is measured at 5100 A (AL 51 n n )- Fo r a review of this method and a comparison 
with RM measurements see IDenney et al.l (j2009l ). where it is shown that for the sufficiently 
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high S/N spectra, the errors associated with sin gle epoch measurements are on the order 
of 0.1 dex (see also IVestergaard fc Peterson! 12006k who estimate errors in Mbh from optical 
and UV scaling relations to be a factor of 4). This method has been increasingly employed 
in studies of AGNs due to the public availability of large samples of high quality spectra 
from surveys like the Sloan Digital Sky Survey (SDSS). Since the slope of the power-law for 
the Rblr— L relation has not been robustly checked at the lowest luminosities, applying the 
single epoch method to large samples of AGN i s an easy way to fin d low mass candidates 
for follow-up studies using RM. The sample of I Greene fc Hoi (120041 ) employs this method 
to identify 19 candidate AGNs from the SDSS with Mbh < 10 6 M . From this sample we 
take SDSS J114008. 71+030711. 4 (hereafter GH08) for a follow-up reverberation mapping 
campaign in an attempt to constrain the properties of AGNs hosting the lowest mass BHs. 
GH08 has a redshift z = 0.0811, which is actually the median of the entire sample. The SDSS 
images show a resolved spiral galaxy with a moderately luminous AGN. The SDSS spectrum 
of GH08 has fairly narrow broad lines typical of narrow line Seyfert 1 (NLS1) type AGN, 
which are typically th ought to hav e lower Mbh but are accreting at close to the Eddington 
limit (for a review see lGrupdl2000l ). While line ratios as measured from the S PSS spectrum 



woul d place GH08 in the HII region of the BPT diagram (see for example iKewley et al. 
there are sev eral detections of an X-ray point source in both ROSAT and Chandra 



( IGreene &: Hoi 120071 ) at the AGN location. HST images presented in iGreene et al.l (120081 ) 
reveal GH08 to have a resolved disk component with spiral arms and an elongated bar. The 
surface brightness as a function of radius is fit with a three component model (disk, bar 
and an AGN) and the AGN luminosity after deconvolution fro m the host is estimat ed to 
be AL 5100 = 5.9 x 10 42 erg s -1 . Using the Rblr— L relation from iBentz et al.l ( )2009al their 
equation 2) and the AGN luminosity estimated from the HST images, we expect GH08 to 
have Rblr ~ 8.7 days. 

In Section 2 we describe the observations, in Sec tion 3 the data analy sis, and in Section 
4 the conclusions. We assume the cosmology used by IGreene fc Hoi (J2004J) (H Q = 100 h = 72 



kms 1 Mpc 1 , Q m = 0.3, Q\ = 0.7) to compare with their results. 



2. Observations &: Data Reduction 



The spectral data for GH08 were taken using the 9.2-m Hobby-Eberly Telescope (HET) 
at McDonald Observatory from December 7, 2009 to June 3, 2010 with a total of 35 obser- 
vations over a period of 178 days. The average time between observations is 6.1 days with a 
median of 5.9. There are 5 instances where the time between useful observations is less than 
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3 days. 

All observations were taken using the Low Resolution Spectrograph (LRS) with the g2 
grism which has a resolving power of R = 1300 and covers a wavelength region from 4300 
— 7300 A. Of these 35 observations, 5 had very low signal to noise (S/N < 2) due to bad 
observing conditions and are not used in the time series analysis. All observations have an 
exposure time of 25 min (providing 30 spectra each with a S/N >3). A fixed slit width of 
2" was used for each observation, but the position angle of the slit varies randomly from 
one observation to the next. A fixed slit position was requested for these queue observations 
but not obtained. The errors associated with a changing position angle will be discussed in 
detail below. 

The spectra were reduced using the standard IRAF routines with the appropriate pa- 
rameters set for the HET LRS instrument. We use an extraction window of 2.5" to maximize 
the S/N while minimizing host contamination from scales larger than about 2 kpc. The spec- 
tra were initially flux calibrated using spectral standard stars taken nightly. On six nights no 
calibration star was observed, usually due to degrading weather conditions, and the standard 
star used for the flux calibration of GH08 on those nights was the one observed closest in 
time to that observation. Wavelength calibration was done using a HgCdAr and Ne spectra 
taken nightly. 

To obtain internal flux calibration between th e spec tra from different nights, the chi- 



square fitting routine of Ivan Groningen &: Wanders! ( 119921 ) was employed. In this algorithm 



the assumption is made that the narrow line flux of the [Olll]A5007 line is constant, since the 
NLR is at a large enough distance that the changing continuum will not effect its line flux 
on these short time scales. Each individual spectrum was then scaled so that the [Olll]A5007 
line in each spectrum had the same integrated flux as the mean spectrum (from averaging 
individual spectra) so that calibrated H/3 measurements could be made. In principle, this 
procedure should have been repeated using the narrow [Sll]AA6716,6731 lines to properly 
calibrate the Ha measurements, but the lines are severely blended and the S/N falls drasti- 
cally just redward of Ha so this was not feasible. Therefore, the results based on the Ha line 
are not as reliable as those for H/3. Other lines originating in the BLR like HelA5876 and 
HellA4868 are not strong enough (or absent) given the S/N to make reliable measurements. 

The rest frame mean spectrum is shown in the top panel of Figure 1 and was ob- 
tained using the highest S/N spectra throughout the run. The most notable features in the 
mean spectrum are the strong Balmer lines, the HelA5876 line, the [Om]A4959,5007 doublet 
and Fell lines flanking H/3. In the mean spectrum we measure the integrated flux of the 
[Om]A5007 line using an extraction region of 80 A (rest frame) about the line centroid and 
find F5007 = 1.67 x 10 _15 erg cm -2 s _1 . The bottom panel of Figure 1 shows the rest frame 
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RMS spectrum of GH08 after the scaling procedure was completed. The [Olll] doublet dis- 
appears as expected and the residual in the H/3 line is the reverberation signal that is to be 
measured. The strong residuals in the Ha, [Nll]AA6548,6584, and [Sn] lines are due to the 
poorer flux calibration in the red part of the spectrum mentioned above. 



3. Data Analysis 



The top panel in Figure 2 shows the light curves for the continuum at 5100 A (rest 
frame). To measure the continuum we take a 140 A wide window centered at 5100 A. This 
window is just redward of the [Om]A5007 line and just blueward of the Fell lines that start 
at ~5200 A (rest frame). The light curves of the integrated line flux for Ha and H/3 are 
shown in the two bottom panels. The measured fluxes for the continuum and Balmer lines 
are listed in Table 1. The light curve statistics are presented in Table 2. The three light 
curves all have a similar global shape, where the variations in the lines are very similar to 
the variations in the continuum. They are however different enough from each other that we 
have confidence that we are not just sampling noise in the seeing conditions or instrument. 
There is roughly a factor of 2.3 overall change in the H/3 and Ha flux, and nearly a factor 
of 2 change for the continuum flux (before subtraction of the constant contribution from 
the host) as shown by R max in Table 2. A possible source of contamination we consider, 
due to the observing procedure used, is that the position angle of the slit was not constant 
throughout the observations. This contamination can manifest as additional narrow line 
flux from extended regions when more of the host galaxy is in the slit, thereby changing the 
relative flux scaling obtained during the inter calibration procedure. Using images of the slit 
position with respect to the host galaxy orientation, we have checked that measured changes 
in flux are not correlated with the slit position angle (thus having more or less of the host 
galaxy in the slit) at different times throughout the observing run. 

For the time series analysis we use two independent methods to cross-correlate the 
continuum and line light curves. The fir st m ethod is the interpol a ted cr oss-correlation func- 
tion (ICCF) of IGaskell fc Sparkd (119861 ) and IGaskell &: Petersonl (119871 ) , as implemented by 



White fc Peterson! ( 



19941 ); see also the review by 



Gaskelj (1994). The second method is the 



is an lm- 



z-transformed discrete correlation function fZDCF) of lAlexander ( 19971) wh ich is 
provement on the discrete correlation function (DCF; Edelson fc Krolikl 1988 ). The ZDCF 
applies Fisher's z transformation to the correlation coefficients, and uses equal population 
bins rather than the equal time bins used in the DCF. In Figure 3 we show the ICCF and 
ZDCF for both the Ha and H/3 lines. The value of the ICCF coefficient for H/3 at the central 
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peak is a modest 0.5 and with a time lag of ~2 days, with an error of +4 and —2 days 
(observed frame). The errors in the time lag determination are estimated using the model- 
independent flux ra ndomization /random subset selection (FR/RSS) Monte Carlo method 
(for more details see lPeterson et al.lll998l ). Our measured value of 2 days is consistent with 
zero time lag within the given errors, meaning that we need finer time sampling to truly 
resolve the reverberation signal, but our measurements constitute an upper limit to the time 
lag. The time lag of zero for Ha could constitute a true lower limit or be spurious due to 
the poor flux calibration as discussed above. 

In Figure 4 we show the velocity profiles of Ha and Hf3 for the mean spectra along with 
the Gaussian line fits used to measure the FWHM. In order to constrain the width of the 
narrow lines we take the [Om]A5007 line as a template and a Gaussian fit gives a FWHM 
of 519 kins" 1 (the effect of instrumental broadening on the line widths and uncertai n ties is 
small). Comparing this to a FWHM of 297 kms -1 as determined by iGreene fc Hoi (120041 ) 
shows that we are only just resolving the narrow lines. In order to constrain the contribution 
of the narrow Balmer lines to the total line profiles we first assume that the [Olll]/H/3 flux 
ratio should be at least 3 b ased on emission line diagnostics for the classification of type 
2 AGN as in iKewley et al.l ( 120061 ). Further, the contribution of the narrow Ha component 
should be roughly 3 times that in the limit of case B recombination. For Ha (top panel) 
we fit 3 narrow line components and one broad. All 3 narrow lines are constrained to have 
the same line widths as [Oin]. The narrow [Nil] lines are forced to have a flux ratio of 2.96 
and the centroids of each line are fixed to the rest frame wavelengths of AA6548,6584. The 
fit falls below the line profile near the peak. This is due to the fact that a simple Gaussian 
function is not quite adequate. The broad component fit should be taller to compensate, 
which would decrease the FWHM. Therefore the FWHM of 737 ± 155 kms -1 we measure for 
Ha can be thought of as an upper limit. The H/3 line (lower panel) is fit with 2 Gaussians, 
a narrow component whose width is modeled after the [Oin] line and one broad component. 
Again the fit falls below the line profile near the peak. We find an upper limit to the 
FWHM of H/3 to be 777 ± 110 kms" 1 . For NLSl-like AGNs a Lorentzian profile can 
pr ovide a better fit to the w ings of the broad component as expected from the analysis 
of iKollatschny fc Zetzll ( 1201 ll ) who show that the broad line shapes can vary from more 
Lorentzian profiles at small FWHM (< 2000 kms" 1 ) to more rectangular profiles at larger 
FWHM (> 6000 kms -1 ). For completeness (but not shown in Figure 4) we also fit both lines 
with a single Lorentzian assuming that the narrow line contribution is negligible. Because 
they are less constrained, the Lorentzians fit better to the overall line profiles. In table 3 we 
show the values for each line as measured from the fits as well as the values measured by 
Greene fc Hoi (120041). For the s ingle Lorentzian fits we find better agreement with original 



measurements of 



Greene fc Hoi (120041 ) . The H/3 line has larger measured FWHM than Ha in 
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all cases. The FWHM of both the H/3 and Ha lines are always larger in the mean sp e ctrum 



than in the RMS spectrum. This is not surprising given the results of ICollin et al.l (120061 ) 
who show that the H/3 line widths are usually 20% broader in the mean spectrum compared 
to the RMS spectrum. We can attribute the overly large widths measured to the fact that 
the lines are just resolved. Therefore, since the true line profiles are not expected to be 
only purely Gaussian or Lorentzian, when estimating the upper limit to Mbh we use an H/3 
FWHM of 703 ± 110 kms -1 which is the average of our two H/3 measurements. 

In order to quantify the contribution of the host galaxy to the AGN luminosity we 
use the analysis by . They used HST/ACS images in two filters (B band — F435W and I 
band — F814W) to decompose the observed intensity distribution of all the images of the 



Greene fc Hoi (120041 ) sample into a combination of components, which includes a central 
unresolved source and some combination of a possible disk, bulge, and a bar. In the case of 
GH08 they do not fit for a bulge, possibly due to the contrast with the strong bar. However, 
the presence of a sm all bulge can not be ruled out from a visual inspection of the images. 



Greene et al.l (120081 ) plot in their Figure 1 the contributions of the AGN, the bar, and the 
disk components as a function of radius, as well as the total surface brightness in the I band 
and the B— I color. We use these components to estimate the AGN contribution to the 
total flux within a circular aperture of radius 2" by integration over the surface brightness 
profiles plotted for the different components of GH08. We find that in the I band the 
AGN contribution is about 20% of the total flux and in the B band it is about 50% of the 



total flux. Using the GH08 AGN luminosity given in iGreene et al.l ( 120081 ) and the HST 
images analysis we described above, we find the galaxy contribution to the luminosity in the 
central 2" of GH08 to be AL 5 i 00 = 7.5 x 10 42 erg s-^Fsioo ~ 1-1 x lO" 16 erg cm" 2 s" 1 A" 1 ). 
Subtracting this contribution from the total luminosity we measure in this study, we estimate 
the AGN luminosity at 5100 A to be about AL 5100 = 1.1 x 10 43 erg s- 1 (F 5100 « 1.6 x 10~ 16 
erg cm -2 s _1 A -1 ). 

We also used the SDSS images of GH08 to estimate the galaxy contribution. We used the 
IRAF isophot package to fit ellipses to the SDSS images of the object and linearly extrapolate 
the radial profile of the elliptical isophotes intensities in the outer 2"— 5" of the image to the 
inner 2" (same as used above with the HST data). The ratio of counts of the measured profile 
inside 2" to the extrapolation from outside 2" is taken as a measure of the AGN fraction 
with respect to the host. We find the AGN to be 46% and 36% of the total luminosity in the 
SDSS g i and r' bands respectively, and we take the average of these two, 41%, as the AGN 
fraction. We note that the AGN fraction may be closer to 46% since the continuum region at 
5100 A is close to the red side of the SDSS g l band. From this extrapolation method we find 
for the AGN that AL 5 i 00 « 7.6 x 10 42 erg s -1 (F 51Q0 « 1.1 x 10~ 16 erg cm" 2 s~ l A^ 1 ) with an 
upper limit of 1.8 x 10 43 erg s -1 (F5ioo = 2.6 x 10~ 16 erg cm -2 s _1 A -1 ). Given the amplitude 
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of continuum variability as shown in Figure 2, these results are in good agree ment with the 



result s we have calculated above from the HST/ACS images analysis done by iGreene et al. 



( 120081 ). In the following we take the AGN luminosity at 5100 A to be 1.1 x 10 43 erg s 1 with 



an an upper limit of 1.8 x 10 erg s . 

Combining all the above measurements and using Equation (1) we find for GH08 that 
M BH = 1.5 x 10 5 M^with an upper l i mit of Mrh 5, 5.8 x 10 5 M ff . The bolometric lu minosity, 
if L bol « 10L 5100 felvis et al.lll994l iMcLure fc Dunlopl hooA iRichards et"aDl2006h . implies 
Lboi/LEdd = 1-5. In Figure 5 we show GH08 in t he Rrt,r— L and Mr h — L planes ( t riangle) , 
with additional sources (circles) comp i led fr om iBentz et al.l ( 120061 ). iBentz et al.l (j2009al ) . 



Bentz et al. 


(2009b) and 


Dennev et al. 


(2010) 



is just within the observed scatter. GH08 deviates more strongly in the M BH — L plot. With 
the upper limits on both L 510 o and M BH GH08 is even more outlying. This is due to the high 
accretion rate and explains the moderate luminosity for such a low Mrh AGN, as indicated 
by the lines of constant Eddington ratio in the bottom panel of Figure 5. 

We can compare the presently determined Mrh with that expected from the M— a* 
relation. This will giv e insight int o the t he coevolution (or lack the reof) of the b u lge in 



AGN with low M BH . iBarth et al.l (120051 ) measure a* for 15 of the IGreene fc Hoi (|2004j ) 
objects (but not GH08) from high resolution spectra and find that most of the objects fit 
well the extrapolation to lower Mrh using the M— a* relation for more massive BHs. As 
these authors mention in their discussion, uncertainties in the analysis rely moderately on 
Mrh determined using the single epoch method and the selection criteria, where the more 
luminous, and hence the more massive BHs are selected. For our a nalysis we loosely assume 
th at (7 [ottt] oc c* (s e e e.g. , iNelson fc Whittldll996l ; iRice et al.ll2006l ). Using the relation given 
by iGiiltekin et al.l ( 120091 ) and a measured ct[oiii] = 126 kms -1 (derived from IGreene fc Ho 



(12004J ) who measure FWHM [0 ni] = 297 kms" 1 ) we would expect M BH ,* « 1.9 x 10 7 M . The 
surprising discrepancy with the measured M BH which is as large as ~1.5 dex probably reflects 
the large uncertainty associated with using the narrow [Olll] line as a proxy for a*, as it has 
been shown that asymmetry in the blue wings and sy stematic blue shifts of the [Olll] lin e 
centroid are observed for many of the NLSl-like objects ( lBorosonll2005l ; lKomossa fc XulbO OT). 
Also, a ex* = 127 kms -1 for GH08 is about twice the average measured a* in the lBarth et al. 



( 120051 ) sample so it is likely that this value is an overestimation but the discrepancy from 



the RM results will still persist. It is still a point of debate as to whether or not AGNs with 



Greene et al. 


2008; 


Mathur et al. 


2011; 


Xiao et al. 


2011) 



This leaves one to speculate on the true impact AGNs with intermediate BHs have in terms 
of feedback on the host properties and the AGNs ability to regulate simultaneous growth of 
a bulge. 



4. Conclusion 



We attempt to measure the central Mbh i n the intermediate mass AGN candidate 
GH08 using the RM technique. GH08 is originally from the sample of iGreene &: Hoi (120041 ) 
and is a low redshift, late type spiral galaxy and hosts a moderate luminosity AGN of the 
NLS1 variety. Our spectroscopic monitoring campaign spans 178 days yielding 30 useful 
observations. From these data we measure the H/3 time lag to be 2 days, with an upper 
limit of 6 days. To further constrain and robustly resolve the time lag will require much 
denser time sampling for this object. Combining the measured r with the H/3 FWHM of 
800 kms -1 , we find M BH = 1.5 x 10 5 M with an upper limit of Mbh < 5.8 x 10 5 M Q . 
This may constitute the lowest Mbh and highest L/LEdd AGN measured from ground based 
observations by RM to date. GH08 deviates from other low mass AGN candidates like NGC 
4395, which notoriously lacks a bulge, in that it has bar and therefore likely has substantial 
vertical structure above the galactic disk (most likely supported near the center by velocity 
dispersion), and therefore is more likely to follow the M— a* relation, which predicts a larger 
BH mass than the one we measure. While relations like M— a* provide evidence for AGN 
feedback and their impact on galaxy evolution, further investigation into low Mbh objects 
like GH08 will constrain deviations from the current paradigm for the coevolution of the 
central SMBH and the host galaxy. 
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Fig. 1.— a. Mean spectrum of GH08; b. RMS spectrum of GH08. 
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Fig. 2. — Light curves for the continuum at 5100 A(top), Ha (middle) and EL3 (bottom). 
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Fig. 3.— Solid line is the ICCF and the dashed line is the ZDCF; Top: Ha CCF; Bottom: 
H/3 CCF. 
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Fig. 4. — Top: The black data points are the Ha mean velocity profile, the dotted lines are 
the Gaussian fit for each component, and the solid line is the total fit to the line profile. 
Bottom: Same as above but for the H/3 mean velocity profile. Both fits fall short at the peak 
due to fitting with pure Gaussians, therefore the measured FWHM can be considered to be 
upper limits. Measured values are presented in Table 3. 
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current best published estimate for the Rrt.r — L relat i on by iBentz et al. f 2009ah (it is not 
a fit to all the shown data p oints - see IDennev et al.l (|2010r) fo r details ); Bot tom: Mrh vs. 
Lfnon- Data points are from Bentz et al. (j2009a ). Bentz et al. ( 2009b ). and Penney et al. 
(120101 ) . The dashed lines are lines of constant Eddington ratio (L bo i/L Edd = 0.01, 0.1, 1.0 



from the top). In both panels the point of GH08 is shown as a red triangle. 
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Table 1. Continuum, Ha and H/3 Light Curves 



Julian Date F ° F Hct F H a 



(1) (2) (3) (4) 



5172 


.99 


1 


.895 


± 





.007 


2 


.668 


± 





.014 





.528 


± 





.005 


5185 


.96 


2 


.214 


± 





.008 


2 


.371 


± 





.013 





.492 


± 





.006 


5190 


.96 


2 


.655 


± 





.008 


2 


.619 


± 





.027 





.498 


± 





.008 


5197 


.93 


2 


.070 


± 





.011 


2 


.187 


± 





.026 





.408 


± 





.013 


5200 


.92 


2 


.618 


± 





.015 


2 


.770 


± 





.034 





.483 


± 





.014 


5208 


.91 


2 


.808 


± 





.012 


2 


.849 


± 





.021 





.595 


± 





.006 


5212 


.89 


2 


.294 


± 





.006 


2 


.671 


± 





.016 





.575 


± 





.004 


5215 


.88 


2 


.354 


± 





.005 


2 


.687 


± 





.019 





.541 


± 





.006 


5221 


.86 


1 


.481 


± 





.005 


1 


.326 


± 





.010 





.253 


± 





.003 


5233 


.95 


1 


.849 


± 





.008 


2 


.304 


± 





.017 





.480 
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.006 


5240 


.82 


2 


.093 
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.007 


2 


.764 


± 





.022 





.490 


± 





.005 
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.80 


2 


.419 
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.006 


2 


.817 


± 





.022 





.506 


± 





.005 
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.79 
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.140 
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.005 


2 


.600 
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.018 





.579 
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.004 
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.005 
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.888 
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.018 





.576 
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.761 
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.858 
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o 


,039 
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,454 


-|- 
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.010 


5259 


.76 


2 


.365 


± 





.008 


2 


.240 


± 





.023 





.444 


± 





.010 
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.88 


2 


.813 


± 





.009 


2 


.876 


± 





.024 





.474 


± 





.007 
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.77 


2 


.293 


± 





.010 


2 


.622 


± 





.022 





.506 


± 





.004 


5272 


.87 


2 
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± 





.006 


2 
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± 





.016 
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± 





.005 
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1 
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2 
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.496 
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1 


.983 


± 
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.454 


± 





.003 


5295 


.66 


1 


.937 


± 





.009 


1 


.744 


± 





.020 





.339 


± 





.007 


5305 


.64 


2 


.126 


± 





.006 


2 


.774 


± 





.019 





.523 


± 





.006 


5315 


.74 


2 


.356 


± 





.012 


2 


.541 


± 





.028 





.471 


± 





.011 


5321 


.72 


2 


.114 


± 





.008 
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.441 


± 





.013 





.526 


± 





.004 
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.72 


2 


.172 


± 





.006 


3 


.020 


± 





.022 





.508 


± 





.005 


5333 


.68 


2 


.241 


± 





.005 


2 


.642 


± 





.015 





.505 


± 





.004 


5341 


.66 


2 


.011 


± 





.010 
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.561 


± 





.029 





.424 


± 





.011 


5350 


.64 


2 


.338 


_|_ 





.005 
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.517 







.014 





.499 


± 





.004 



Note. — Col.(l): J.D. -2450000 days; Col.(2): Continuum 
flux (10~ 16 erg cm" 2 s" 1 A" 1 ); Col.(3) & (4): Line flux (10~ 14 
erg cm~ 2 s -1 ). The length of each exposure is 0.0174 days. 



-20 - 



Table 2. Continuum, Ha and H/3 Light Curve Statistics 



Cont./Linc 


Fniin 


Fmax 


Rmax 


Favg 


<7F 


Fvar 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Cont. at 5100 A 


1.48 


2.81 


1.90 


2.23 


0.298 


0.134 


Cont. at 1500 A 


0.38 


1.71 


4.50 


1.13 


0.298 


0.287 


after host subtracted 














H/3 


0.25 


0.59 


2.35 


0.48 


0.069 


0.141 


Ha 


1.33 


3.02 


2.28 


2.53 


0.359 


0.142 



Note. 

Line flux in units of 10 -14 erg cm~ 2 



erg cm 2 s 1 



A 



-l. 



Continuum flux in units of 10 16 
is the standard deviation; F var is the fractional variation as defined by 



Rodriguez-Pascual et al.l ([19971 ). 
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Table 3. [Om], Ha and H(3 Measured Velocities 



Line 


Gaussian FWHM 


Lorentzian FWHM 


G&H04 




(km s _1 ) 


(km s- 1 ) 


(km s _1 ) 


(1) 


(2) 


(3) 


(4) 


[Om] 


519 ± 11 


454 ± 19 


297 


Ha 


737 ± 155 


540 ± 6 


591 


H/3 


777 ± 111 


629 ± 12 





Note. — Line widths for [Om], Ha and H/3 in km s _1 . The 
Gaussian fits include both a narrow and broad component for 
the Balmer lines and are shown in Figure 4. The Lorentzian fits 
are for a single component where the narrow line contribution 
is as sumed to be negligibl e. The values in the final column are 
from iGreene fc Hoi d2004h . 



